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^ Abstract 

)^ ' 

, ^ We derive the expression for the residual interaction to be used in the framework of the Second RPA with density-dependent forces. 

^> The adopted procedure is based on a variational approach. It is found that the residual interaction to be used in RPA and beyond 
RPA matrix elements is not the same as far as the rearrangement terms are concerned. A detailed derivation of the matrix elements 
coupling 1 particle- 1 hole with 2 particle-2 hole and 2 particle-2 hole among themselves has been done and the corresponding 
rearrangement terms are shown here. This formal result indicates that both the currently used prescriptions, namely (i) using the 
same type of rearrangement terms in RPA and beyond RPA matrix elements or (ii) neglecting the rearrangement terms in beyond 

'"pi RPA matrix elements, are not correct. 

^Keywords: Variational method, second random-phase approximation, residual interaction 
r3 PACS: 21.30.-x, 21.60.-n, 21.60.Jz 
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The RPA is a widely used microscopic approach adopted to 
study collective modes which can be interpreted in terms of vi- 
brations. In nuclear physics, the employment of this model with 
effective density-dependent interactions allows to successfully 
describe giant resonances and low-lying excitation modes. The 
second RPA (SRPA) is a natural extension of the RPA where 
2 particle-2 hole (2p-2h) excitations are included together with 
the usual RPA 1 particle- 1 hole (Ip-lh) configurations, provid- 
ing in this way a richer description of the excitation modes. On 
the other hand, the manifestation of new exotic modes in unsta- 
ble nuclei has produced an increasing interest in the context of 
nuclear structure and the inclusion of complex configurations 
coupled among themselves may play a crucial role in the de- 
scription of these phenomena. Moreover, some low-energy ex- 
citations like the first 0^ and 2^ states in magic nuclei are not 
well reproduced by the standard RPA model (they are predicted 
too high in energy) because beyond Ip-lli configurations are 
needed to describe them. Finally, SRPA would allow to study 
in a proper way the double phonon excitation modes that are 
experimentally well known in nuclei. In all these cases, inter- 
esting information may be provided by the application of SRPA. 

The SRPA equations are well known since several decades 
(see, e.g.. 111]). However, up to very recently, they have never 
been fully and self-consistently solved due to the heavy numer- 
ical effort that they require. Some approximations have been 
adopted in the past, namely the SRPA equations have been re- 
duced to a simpler second Tamm-Dancoff model (see for in- 



stance 10 SB Si) and/or the equations have been solved with 
uncorrected 2p-2h states, i.e. the residual interaction terms in 
the matrix that couples 2p2h configurations among themselves 
have been neglected (diagonal approximation) 
lll[l2|l- Very recently, the SRPA equations have been solved for 



some closed-shell nuclei using an interaction derived from the 
Argonne VI 8 potential (with the Unitary Correlation Operator 
Method) 11 1311 . for small metallic clusters in the jellium approx- 
imation lll4|] and for the magic nucleus '^O with the Skyrme 
interaction OlSll . 

The most currently employed phenomenological interac- 
tions in mean field approaches are density-dependent forces of 
Skyrme or Gogny type. It is well known that, with density- 
dependent forces, the residual interaction used to evaluate the 
RPA matrices A and B contains rearrangement terms gener- 
ated by the derivative with respect to the density of the mean 
field hamiltonian (second derivative of the energy density func- 
tional). Moreover, rearrangement terms are present also in the 
single particle energies appearing in the diagonal part of the 
matrix A. Due to the density dependence of the interaction, the 
mean field is changed and the single particle energies are thus 
also modified by the variations of the density. 

When dealing with the SRPA problem with density- 
dependent forces, a first formal aspect to consider is the de- 
termination of the residual interaction which has to be used in 
beyond RPA matrix elements. To our knowledge, this aspect 
has not yet been clarified in the literature. 
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In Ref. lfl6ll the issue of the rearrangement terms has been 
studied in the context of shell model calculations while in Ref. 
UtIi a prescription is introduced for the rearrangement terms 
appearing in matrix elements beyond the standard RPA ones. 
However, in more recent calculations lIlSll . the same authors 
have not actually used that prescription and have not included 
those rearrangement terms. 

This work differs from that of Ref. UtIi for the following 
reasons. The ground state of the system is described as a su- 
perposition of Ip-lh and 2p-2h configurations build on top of 
the Hartree-Fock (HF) state, (see Eq. (|2|). The coefficients de- 
scribing these configurations are used as variational parameters 
and the SRPA equations are obtained, as shown in Ref. I.l9i1 . by 
minimizing the expectation value of the Hamiltonian with re- 
spect to them. On the contrary, in Ref. [ 17] the 2p-2h degrees 
of freedom are not taken into account explicitly. In the case of 
density dependent forces, the RR\ equations can be derived in 
the small amplitude limit of the Time Dependend Hartre-Fock 
(TDHF) method which allows to define the residual interaction 
as the second derivative of the energy density functional with 
respect to the density. This implies that the rearrangement terms 
appear and they are related to the first and second derivative of 
the interaction with respect to the density. Therefore, differently 
from Ref. fT?], we consider an expansion of the interaction up 
to the second order in the density. This allows obtaining the 
correct definition of the rearrangement terms in the RPA ma- 
trices. Finally, the analysis done in (vf\ deals only with the 
rearrangement terms in the matrix elements between the \p-\h 
and 2p-2h configurations while the coupling of the 2p-2h con- 
figurations among themselves is not studied. 

In this work, we discuss the issue of the residual interaction 
in SRPA using the variational derivation of the SRPA equations 
proposed by Providencia lfl9ll . extending this derivation to the 
case of a density dependent force. As done in 1.19.1 . we express 
the ground state \V > as 



The one body density p matrix can be expanded around the 
HF density p*"' 



p^p -< *I'|fl^fl„|>I' >=< ^\e^' a'^aae^\<^ > 
=< 0|(1 H- 5"' + + ...)a^Ml +S + ^S^ + ...)|<1) > 



< ^ajjUaS ^ + S ''a^OaS + '^aj^a„\<i> >= p^^ -n 6pa/s (6) 

where we have truncated at the quadratic terms. Hence, the 
variation of the density is given by the sum of a linear 5p*'^ and 
a quadratic 6p^-^ contributions written as; 



and 



expS\<^> , 



(1) 



In particular, one can find that: 

Mil - = ; sp^'J, = c,, ■ sp['l = c;, ■ 



(7) 



(8) 



(9) 



(10) 



(11) 



(12) 



where > is the HF ground state and 5 is a linear superposi- 
tion of Ip-lh and 2p-2h configurations 



ph pp'hh' 



2 ^ --pp'hh'^p^jj'Clh'Clh 
pp'hh' 



(2) 



considering thus a different and more general wave function 
than the one used in [17]. 

Let us denote with h, i, j, k and I hole states and with m, n, p 
and q particle states. We introduce: 



(3) 



We consider the case of a density-dependent interaction and 
the Hamiltonian is thus written in the following form: 



ap afiyS 



where: 

Val3yd(p) = Val3ys(p) - Val36y(p) 



(4) 



(5) 



The mean value of the Hamiltoninan in >, < H >, can be 
written as in Eq. (8) of Ref. HI]: 

< H >=< 0|//|<D > + J](C;,.i™(p) + C„uAi,„(p)) + 

mi 

X (C:n„uV'nnij(p) + C,nnijVij>nn(p)) + F'^^^ , (13) 
/'< j,m<n 

where F*^^ is the sum of several quadratic contributions in C 
and C (see Eq. (8) of Ref. IH] for details). 
In the above equation: 



k 



(14) 



which, in cases where the interaction is not density-dependent, 
define the single particle energies. The A matrices appearing in 
the SRPA equations are defined by 



^mi,pk — 



6^ <H> 



SCIJC 



pk 



(15) 



^mi,pqkl ~ 



6"^ <H> 



5C*:5C 



pqkl 



= A 



^mnij,pqkl ~ 



6^ <H> 



12, 



= A 



22, 



(16) 



(17) 



which are evaluated at C = C* = C = C* = 0. The corre- 
sponding B matrices are obtained by substituting in the previ- 
ous equations the derivative with respect to C and C with the 
derivative with respect to C* and C*, respectively. For example, 
the RPA B matrix is 



B 



mi.pk ~ 



6^ <H> 



SC*:6C 



pk 



(18) 



When the interaction is not density-dependent, only the F*^' 
term of Eq. ( fTlT l gives some contributions, since the others 
are not quadratic in the C's coefficients. The standard SRPA 
matrices are thus obtained (see Eqs (9)-(12) of Ref. LiSJ). If, 
on the contrary, the interaction is density-dependent, rearrange- 
ment terms appear and they are given by the first three terms of 
Eq. ( fTST l. In the following we focus our attention only on the 
rearrangement terms that have to be included in SRR\ matrices. 

As mentioned above, in the case of density dependent forces, 
RPA equations can be obtained in the small amplitude limit 
of TDHF and the single-particle energies contain only the first 
derivative of the interaction with respect to the density, while 
both the first and the second derivatives appear in the residual 
interaction. We thus expand V around the HE density p^^^ and 
keep up to quadratic terms : 



ab 



4z 

abed 



5V» 



./SyS 



Spah 



Spa 



'PyS 



SpahSpcd 



SpahSpcd ; 



where 



(19) 



(20) 



It can be shown that, by using the above definition for the 
RPA matrices and the previous expansion for the interaction, 
the usual RPA rearrangement terms are found. 

Eor example, for the An matrix defined in Eq dTSl ) we find 



^mi,pk — ^ki^mp ^mp^ki "t" ^mkip, 

with the single particle energies given by 



^ kt 



SVkk'kk' 



Spba 



Pha, 



(21) 



(22) 



p=p"' 



where pba - 4'1,4'a and if/ are the HE single particle wave func- 
tions. The rearrangement contribution in Eq. ( l22l l comes out 
from the derivative of the first term appearing in the right side 



of Eq. ( fTSl i; in the derivation of Eq. (l22l) we have used that 
Phh' ~ ^i''^'- ^1' the residual interaction is 



'^mkip 
SVhkhp 



Vrnkipip'^"^) + 2 



6V, 



mhih 



6pkp 



Pkp 



p=p'" 



Spn 



P, 



hh' 



hh'hh' 



SpmiSpkp 



Pmipkp, 

(23) 



where the first term is the original interaction while all the oth- 
ers are the rearrangement terms. The first two come out from 
the A quantities appearing in Eq. ( fT3T l, while the last one derives 
from the mean value of the Hamiltonian in the HE state, i.e. the 
first term of Eq. ( fT3T l. 

Let us now focus our attention on the rearrangement terms 
appearing in the residual interaction of the matrix elements be- 
yond the usual RPA. Eirst of all, it comes out that no rearrange- 
ment terms appear in the matrix A22, Eq. (fTTl l. except for the 
contributions 



^(rearr) 
mnij.pqkl 



1 

- +-5, 



qSikdjl ^ 
kk' 

--SmpSikSjl 
^ kk' 



SVkk'kk' 



Sppm 

SVkk'kt 



1 



^ kk' 
1 ^ 

^ kk' 



5pqn 
SVkk'kk' 



Ppm 



Pqn 



Spki 
SVkk'kk' 



Spij 



Jp=pl" 



Pki 



Plj 



(24) 



generated by the first term of Eq. (fT3T l. The above terms, to- 
gether with the A quantities, give the correct single particle en- 
ergies in the case of a density-dependent interaction, being thus 
completely consistent with the single particle energies appear- 
ing in the RPA matrix (1211 1. No rearrangement terms appear in 
the B22 matrix. 

Let us consider now the matrices coupling 1 p-lh with 2p-2h 
configurations: the same kind of contributions discussed just 
above is found, that together with the first term appearing in Eq. 
10 of ref. 1,19.] gives a vanishing term due to the HE condition 



ekp = ^*p(p<°') + ^ 

kk' 



kk' 



Sppk 



Ppk 



0. 



(25) 



In addition to that, we have only one type of rearrangement 
terms in the residual interaction, given by: 



^(rearr) 
mi, pqkl 



6Vkln 



6ph 



Pin. 



A similar term is found in the B\2, namely 



(rearr) 
mi, pqkl 



6Pm 



(26) 



(27) 



We recall that, in the case of a two-body non density- 
dependent interaction, the B22 and B\2 are found to be zero. 
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On the contrary, in the present case Bn has a non vanishing 
contribution coming from the rearrangement terms. This can 
be easily understood since the density-dependent part of the in- 
teraction mimics a three-body force and takes thus into account 
more comphcated correlations. 

Some conclusions can be thus drawn for the residual interac- 
tion of the matrix elements beyond RPA. A first kind of rear- 
rangement terms appears in the matrix A22 giving the right def- 
inition of the single particle energies, consistently with RPA. 
The same kind of terms appears also in the A 12 matrix which 
finally gives a vanishing contribution due to the HF condition 
Eq. ( l25T l. No rearrangement terms appear instead in the residual 
interaction of the A22 and B22 matrices. The B22 matrix is thus 
still equal to zero. A comment can be added concerning this 
last result, that, as a matter of fact, introduces an asymmetry 
between the An and A22 matrices, as far the residual interac- 
tion is concerned. This can be traced back to the fact that the 
interaction depends only on the one-body density. A different 
situation would occur in the case of more general forces, for 
example depending also on the two-body density matrix. 

In the matrices A 12 and B12, rearrangement terms of the resid- 
ual interaction are found and they are given by Eqs. ( |26] | and 
{Tl\ . respectively. Their expression is different from the one of 
the RPA matrices. Therefore, the same definition of the residual 
interaction as in standard RPA cannot be adopted in the matrix 
elements describing the coupUng between 1 particle- 1 hole with 
2 particle-2 hole and 2 particle-2 hole among themselves. At 
the same time, completely neglecting the rearrangement terms 
in the matrix elements beyond RPA is not correct. A quantita- 
tive study of the terms ( |26] | and (l27T i could be useful and it will 
be done in a future work. 
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